Measuring Statistical Isotropy of CMB 

Anisotropy 



Tarun Souradeep\ Amir Hajian^'^ and Soumen Basak^ 

^ Inter- University Centre for Astronomy and Astrophysics (lUCAA), 
Post Bag 4, Ganeshkhind, Pune 4^1 007, India. 
E-mail: tarun@iucaa.ernet.in;soumen@iucaa.ernet.in 
^ Department of Physics, Jadwin Hall, Princeton University, Princeton, NJ 08542. 
^ Department of Astrophysical Sciences, Peyton Hall, Princeton 
University, Princeton, NJ 08544- 
E-mail: ahajian@princeton. edu 



Abstract 

The statistical expectation values of the temperature fluctuations and polariza- 
tion of cosmic microwave background (CMB) are assumed to be preserved under 
rotations of the sky. We investigate the statistical isotropy (SI) of the CMB maps 
recently measured by the Wilkinson Microwave Anisotropy Probe (WMAP) using 
the bipolar spherical harmonic formalism proposed in Hajian & Souradeep 2003 for 
CMB temperature anisotropy and extended to CMB polarization in Basak, Hajian 
& Souradeep 2006. The Bipolar Power Spectrum (BiPS) had been measured for the 
full sky CMB anisotropy maps of the first year WMAP data and now for the recently 
released three years of WMAP data. We also introduce and measure directional sen- 
sitive reduced Bipolar coefficients on the three year WMAP ILC map. Consistent 
with our published results from first year WMAP data we have no evidence for 
violation of statistical isotropy on large angular scales. Preliminary analysis of the 
recently released first WMAP polarization maps, however, indicate significant viola- 
tion of SI even when the foreground contaminated regions are masked out. Further 
work is required to confirm a possible cosmic origin and rule out the (more likely) 
origin in observational artifact such as foreground residuals at high galactic latitude. 
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1 Introduction 



In standard cosmology, CMB anisotropy signal is expected to be statistically 
isotropic, i.e., statistical expectation values of the temperature fluctuations 
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AT{q) are preserved under rotations of the sky. In particular, the angular 
correlation function C{q, q') = {AT{q)AT{q')) is rotationally invariant for 
Gaussian fields. In spherical harmonic space, where AT{q) = J2im. '^imXimiQ) 
the condition of statistical isotropy (SI) translates to a diagonal {aimO*ii^,) = 
Ci6ii'6mm' where C/, is the widely used angular power spectrum of the CMB 
anisotropy. Statistical isotropy of the CMB sky is essential for Ci to be a 
complete description of (Gaussian) CMB anisotropy and, hence, an adequate 
measure for carrying out cosmological parameter estimation of 'standard' (SI) 
model. Hence, it is crucial to be able to determine from the observed CMB sky 
whether it is a realization of a statistically isotropic process, or not. The detec- 
tion of statistical isotropy (SI) violations in the CMB signal can have exciting 
and far-reaching implications for cosmology and the cosmological principle. 
For example, a generic consequence of cosmic topology is the breaking of sta- 
tistical isotropy in characteristic patterns determined by the photon geodesic 
structure of the manifold as probed by the CMB photons traveling to us from 
the surface of last scattering over a distance comparable to the cosmic hori- 
zon, Rh dHH;!!)- Mildly anisotropic cosmological models predict charateristic 
patterns hidden in the CMB sky. On the other hand, SI violation could also 
arise from foreground contamination, or, artifacts of observational and analysis 
techniques. 



The CMB measurements of the Wilkinson Microwave Anisotropy Probe ( WMAP) 
are consistent with predictions of the concordance ACDM model with (nearly) 
scale-invariant and adiabatic fluctuations expected to have been generated 
during an inflationary epoch in the early universe (0; IE S S S IE Hi! ). 
After the first year of WMAP data, the SI of the CMB anisotropy [i.e. rota- 
tional invariance of n-point correlations) has attracted considerable attention. 
Tantalizing evidence of SI breakdown (albeit, in very different guises) has 
mounted in the WMAP first year sky maps, using a variety of different statis- 
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The three- year WMAP maps are consistent with the first-year 
maps up to a small quadrupole difference. The two additional years of data 
and the improvements in analysis has not significantly altered the low multi- 
pole structures in the maps (|lOl). Hence, 'anomalies' are expected to persist 
at the same modest level of significance but are now less likely to be arti- 
facts of noise, uncorrected systematics, or the analysis of the first year data. 
The cosmic significance of these 'anomalies', however, remains debatable be- 
cause of the aposteriori statistics employed to ferret them out of the data. 
More importantly, what is missing is a common, well defined, mathematical 
language to quantify SI (as distinct from non Gaussianity) and the ability to 
ascribe statistical significance to the anomalies unambiguously. We employ a 
well defined, mathematical language of Bipolar harmonic decomposition of the 
underlying correlation to quantify SI that can ascribe statistical significance 
to the anomalies unambiguously. 
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The observed CMB sky is a single realization of the underlying correlation, 
hence the detection of SI violation or correlation patterns pose a great observa- 
tional challenge. For statistically isotropic CMB sky, the correlation function 

C{ni,h2) = C{hi-n2) = dnC{nni, 7^n2), (1) 

where TZn denotes the direction obtained under the action of a rotation TZ on 
n, and dTZ is a volume element of the three-dimensional rotation group. The 
invariance of the underlying statistics under rotation allows the estimation of 
C{n- h') using the average of the temperature product AT{h)AT{n') between 
all pairs of pixels with the angular separation 9. In the absence of statistical 
isotropy, C(n, h') is estimated by a single product AT(n) AT(n') and hence 
is poorly determined from a single realization. Although it is not possible to 
estimate each element of the full correlation function C(n, n'), some measures 
of statistical anisotropy of the CMB map can be estimated through suitably 
weighted angular averages of AT{n)AT{n'). The angular averaging procedure 
should be such that the measure involves averaging over sufficient number 
of independent 'measurements', but should ensure that the averaging does 
not erase all the signature of statistical anisotropy. We proposed the Bipolar 
Power spectrum (BiPS) ki {i = 1,2,3, . . .) of the CMB map as a statistical 
tool of detecting and measuring departure from SI (B Q and reviewed in 
this article in sec. 2. The non-zero value of the BiPS spectrum imply the break 
down of statistical isotropy 

STATISTICAL ISOTROPY =^ = y£ ^ 0. (2) 



The BiPS is sensitive to structures and patterns in the underlying total two- 
point correlation function (|35l : bfii) . The BiPS is particularly sensitive to real 
space correlation patterns (preferred directions, etc.) on characteristic angular 
scales. In harmonic space, the BiPS at multipole i sums power in off-diagonal 
elements of the covariance matrix, {aimdi'm'), in the same way that the 'angular 
momentum' addition of states Im, I'm' have non-zero overlap with a state with 
angular momentum |/ — /'| < i < I + I'. Signatures, like aim and ai+nm being 
correlated over a significant range / are ideal targets for BiPS. These are typical 
of SI violation due to cosmic topology and the predicted BiPS in these models 
have a strong spectral signature in the bipolar multipole i space ^3l\\. The 
orientation independence of BiPS is an advantage for constraining patterns 
(preferred directions) with unspecified orientation in the CMB sky such as 
that arising due to cosmic topology or, anisotropic cosmology (jish. 



The results of WMAP are a milestone in CMB anisotropy measurements since 
it combines high angular resolution, high sensitivity, with 'full' sky coverage al- 
lowed by a space mission. The frequency coverage allows for WMAP CMB sky 
maps to be foreground cleaned up to / ~ 100 jT^H^. The CMB anisotropy 
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map based on the WMAP data are ideal for testing statistical isotropy. Mea- 
surement of the BiPS on CMB anisotropy maps based on the first year WMAP 
data are consistent with statistical isotropy (4^ 4l[ li^. BIPS analysis on the 



recent WMAP-3yr data are consistent with SI and further, indicate that BiPS 
of the three years maps show an improvement in SI - the deviations are smaller 
and fewer (|43l). 



First 'full-sky' CMB polarization maps have been recently delivered by WMAP flllh 
The statistical isotropy of the CMB polarization maps are an independent 
probe of statistical isotropy. Since CMB polarization is generated on at the 
surface of last scattering, violations of statistical isotropy are pristine cosmic 
signatures and more difficult to attribute to the local universe. The bipolar 
power spectrum has been defined and implemented for CMB polarization and 



show great promise (j4J). 



2 The Bipolar representation of SI violation 



Although it is not possible to estimate each element of the full correlation 
function C{hi, 11,2), some measures of statistical isotropy of the CMB map can 
be estimated through suitably weighted angular averages of AT(?7.i)AT(n2). 
The angular averaging procedure should be such that the measure involves 
averaging over sufficient number of independent measurements to reduce the 
cosmic variance, but should ensure that the averaging does not erase all the 
signature of statistical anisotropy. The Bipolar power spectrum (BiPS) is a 
measure of statistical isotropy proposed by us in ref. |35i) 

= {21 + If j dQn, J dQn, J dUx^n) C{nn,, TZn^)]' (3) 



that also has another important desirable property of being independent of the 
overall orientation of the sky. In the above expression, C{TZni, TZn2) is the two 
point correlation at TZhi and TZn2 which are the coordinates of the two pixels 
hi and 712 after rotating the coordinate system by element TZ of the rotation 
group, is the trace of the finite rotation matrix in the ^M-representation 
X^(7^) = Z)m=-£ -Omm('^) which is called the characteristic function, or the 
character of the irreducible representation of rank I. For a statistically isotropic 
model, C(ni, 112) is invariant under rotation, and therefore C{TZhi, 7lh2) = 
C{ni, 722) and then the orthonormality of x^{^)j implies the condition for SI, 
= KP6eo stated in eq.(2). 

The harmonic space representation BiPS provides complementary understand- 
ing and allows for rapid numerical computation. Two point correlations of 
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CMB anisotropy, C(ni, 71,2), are functions on S'^ x S'^, and hence can be ex- 
panded as 

C(ni,n2)= ^ ^^m|M2^/m (^1' ^2). (4) 

Here A^f^^^^i^ are the Bipolar Spherical harmonic (BipoSH) coefficients of the 

expansion and Y^'|j^(ni, 77,2) are bipolar spherical harmonics. Bipolar spherical 
harmonics form an orthonormal basis on 5*^ x S"^ and transform in the same 
manner as the spherical harmonic function with ^, M with respect to rotations 
(0). Consequently, inverse-transform of C(ni, 112) in eq. (4) to obtain the 
BipoSH coefficients of expansion is unambiguous. 

Most importantly, the Bipolar Spherical Harmonic (BipoSH) coefficients, ^£jvf|Zi 
are linear combinations of off-diagonal elements of the harmonic space covari- 
ance matrix, 

mim2 

where C^^^i^^^ are Clebsch-Gordan coefficients. This clearly shows that A^f^^i^^^ 
completely represent the information of the covariance matrix. 

Statistical isotropy implies that the covariance matrix is diagonal, {aimof,,^,) = 
Ci Sii'Smm' and hence the angular power spectra carry all information of the 
field. When statistical isotropy holds BipoSH coefficients, A^j^^n,, are zero 
except those with i = 0, M = which are equal to the angular power spectra 
up to a (— 1)'(2/ + 1)^/^ factor. Therefore to test a CMB map for statistical 
isotropy, one should compute the BipoSH coefficients for the maps and look 
for nonzero BipoSH coefficients. Statistically significant deviations of BipoSH 
coefficient of map from zero would establish violation of statistical isotropy. 

It is impossible to measure all Ai^^ individually from the single CMB sky map 
because of cosmic variance since they form an equivalent representation of a 
general two point correlation function. There are several ways of combining 
BipoSH coefficients that serve to highlight different aspects of SI violations. 
Further, one can exploit the freedom to smooth the map in real space by ap- 
propriate symmetric kernels W{h-h') to target specific regions of the multipole 
space by the corresponding window function Wi. 

2.1 Bipolar Power spectrum (BiPS) 

Among the several possible combinations of BipoSH coefficients, the Bipolar 
Power Spectrum (BiPS) has proved to be a useful tool with interesting features 
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f)35f l. BiPS of CMB anisotropy is defined as a convenient contraction of the 
BipoSH coefficients 
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where Wi is the window function that corresponds to smoothing the map in 
real space by symmetric kernel to allow targeting specific regions of the mul- 
tipole. BiPS is interesting because it is orientation independent, i.e. invariant 
under rotations of the sky. SI condition implies a null BiPS, i.e. = for 
every ^ > 0, {ki = Kodm). Non-zero components of BiPS imply break down 
of statistical isotropy, and this introduces BiPS as a measure of statistical 
isotropy. It is worth noting that although BiPS is quartic in aim, it is designed 
to detect SI violation and not non-Gaussianity 
estimator of BiPS is given by 



4l|;|4a). An un-biased 
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where is the bias that arises from the SI part of the map. The bias is given 
by the angular power spectrum, Ci that can be estimated from the map (done 
in the results quoted here), or, correspon ding to the best fit theoretical model 



(done for WMAP-1 analysis) (|23 |4l|; 143; |4 



2.2 Reduced Bipolar Coefficients 



The BipoSH coefficients can be summed over / and /' to reduce the cosmic 
variance, 

oo (.+1 

= E E ^iM\ii'- (8) 

1=0 V=\t-l\ 



In any given CMB anisotropy map, A^m would fluctuate about zero. A se- 
vere breakdown of statistical isotropy will result in huge deviations from zero. 
Reduced bipolar coefficients are not rotationally invariant, hence they assign 
direction to the correlation patterns of a map. An interesting way of visualizing 
these coefficients is to make a Bipolar map from Aim 

oo I 

6>(^)=E E A^MYm{n). (9) 

£=0 M=-t 
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Fig. 1. Top: A bipolar map generated from bipolar coefficients, Aim, of 3-year 
ILC map. Middle: bipolar map based on 1-year ILC map and Bottom: differ- 
ences between the two maps (note the scales). The top map (ILC-3) has smaller 
fluctuations comparing to the middle one (ILC-1) except for the hot spot near 
the equator. Differences between these two maps mostly arise from a band around 
the equator in bipolar space. Both ILC maps are smoothed by a band pass filter, 
W^{lt = 2,ls = l0). 

The symmetry Aim = (~1)^^^|-m of reduced bipolar coefficients guarantees 
reality of B(n). The bipolar map based on bipolar coefficients of ILC-3 is 
shown on the top panel of Fig. 1. The map has small fluctuations except 
for a pair of hot and cold spots near the equator. To compare, we have also 
made a bipolar map of 1-year ILC map (ILC-1) from bipolar coefficients of 
ILC-1 (middle panel of Fig. 1). The difference map (Fig. 1 (bottom)) shows 
that differences between these two maps mostly arise from a band around the 
equator in bipolar space. As it is seen in Fig. 1, the bipolar map of ILC-3 
has less fluctuations comparing to that of ILC-1. This is because almost all of 
A^m's of ILC-3 are smaller than those of ILC-1 {i.e. are closer to zero). 

Figure 2 compares A^m and Kg of three year WMAP ILC map against an 
average of 1000 simulations of statistically isotropic maps with same angular 
power spectrum. As it can be seen many spikes presented in Aim's of ILC-1 
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Fig. 2. Left: Real part of ^£m's of ILC-3 (red square points) and ILC-1 (blue stars) 
for a VF;'^(10,2) filter that roughly keeps multipoles between 2 and 15. £ and M 
indices are combined to a single index n = i{i + 1) + M + 1. The blue dotted lines 
define 1-a error bars derived from 1000 simulations of SI CMB anisotropy maps. 
Almost all of A^^f 's of ILC-3 are smaller than those of ILC-1 which means ILC-3 
is more consistent with statistical isotropy. Right: The BiPS of ILC-3 (red square 
points) and ILC-1 (blue stars) for a iy^'^(10,2) filter is consistent with statistical 
isotropy The green dotted lines define 1-a error bars derived from 1000 simulations 
of SI CMB anisotropy maps. 

have either disappeared or reduced in ILC-3 {e.g. those around n = 20, 40 and 
a big spike at n = 111). More detailed description of analysis and results are 
given in ref. JiT 



The rotational properties of Aim suggests defining another rotationally in- 
variant bipolar power spectrum in Di = J^m AimA}^^ - the same way aim 
are combined to construct the angular power spectrum, C;. The measurement 
of reduced Bipolar Power Spectrum (rBiPS), Di, on CMB maps is ongoing. 
We also defer a comparison of the two orientation insensitive measures of SI 
violation Kg (BiPS) and Di (rBiPS). 



3 BiPS of CMB polarization maps 



One of the firm predictions of this working 'standard' cosmological model is 
linear polarization pattern {Q and U Stokes parameters) imprinted on the 
CMB at the last scattering surface. A net pattern of linear polarization is 
retained due to local quadrupole intensity anisotropy of the CMB radiation 
Thomson scattering on the electrons at at the last scattering surface. 

Recently WMAP has provided the first "full" sky CMB polarization maps. 
The wealth of information in the CMB polarization field will enable us to 
determine the cosmological parameters and test and characterize the initial 
perturbations and inflationary mechanisms with great precision. Cosmological 
polarized microwave radiation in a simply connected universe is expected to 
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BiPR for WUAPSrd (Reduced Foreground) 




Fig. 3. BiPS for polarization maps(nside = 512) with polarization mask using 
band-pass filter Wi^{ls = 10, = 5). Red dots show the BiPS for WMAPSrd fore- 
ground cleaned polarization map with polarization mask after bias subtraction, 
green dots show the BiPS for 100 SI maps with polarization mask after bias sub- 
traction and Blue lines show the 1 — cr of the cosmic variance of BiPS for 100 SI 
maps with polarization mask. 

be statistically isotropic. This is a very important feature which allows us to 
fully describe the field by its power spectrum that can have profound theo- 
retical implications for cosmology. Statistical isotropy (SI) can now be tested 
with CMB polarization maps over large sky fraction. Importance of having 
statistical tests of departures from SI for CMB polarization maps lies not only 
in interesting theoretical motivations but also in testing the cleaned CMB 
polarization maps for residuals from polarized foreground emission. 

The coordinate-free description of CMB polarization decomposes the two 
kinds of polarization pattern on the sky based on their different parities. In 
the spinor approach, the even parity pattern is called the i?-mode and the 
odd parity pattern the i?-mode. Hence the CMB sky maps are characterized 
by a triplet of random scalar fields: X{n) = {T(n), E{h), B{n)}. Statistical 
properties of each of these fields can be characterized by A^-point correlation 
functions, {X{hi)X{h2) ■ ■ ■ X{hn)) ^ ■ 

Gaussian CMB sky is completely described by two-point correlation functions 
of X{h), or equivalently, the corresponding spherical harmonic coefficients 

C^^'(n,n') = {X{n)X'{n')y{ataif::,). (10) 



^ For cut-sky, E{n) and B{n) mode decomposition is not unique (|46l : 47). But since 
mixing is linear there always exist two linearly independent modes. It is possible to 
formulate the SI of these linear independent modes. 
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It is possible to define generalized Bipolar Spherical Harmonic (BipoSH) co- 
efficients, Afj^f^i_^i^, that are linear combinations of off-diagonal elements of the 
{a^^aff^,) covariance matrix, 

aXX' _ / X *X' \(_i\ni2MM /..N 

mim2 

where Cf^_^i^^^ are Clebsch-Gordan coefficients. Incomplete sky coverage in- 
duces a contamination of E-mode of polarization by its B-mode and vice- versa. 
Then the modified temperature and polarization fields is related to their ac- 
tual values of full sky coverage by a window matrix ( 
are basically window functions for temperature and polarization in harmonic 
space. It can be shown that the estimated BipoSH coefficients are in fact linear 
combinations of that for full-sky CMB maps, 

^eM\ii' — X! ^e^nih ^eMihh ■ (l^) 
e'M'hh 

Here bold-faced A^j^^^, and A^^^i^^^^ are the column matrices corresponding to 
estimated and true IBipoSH coefficients respectively, for the auto and cross- 
correlations {TT, TE, TB, ET, EE, EB, BT, BE, BB) of temperature anisotropy 
and polarization. The elements of the matrix N^/M'«ii2 depend on Clebsch- 
Gordan coefficients and window functions in harmonic space. Hence, the true 
BipoSH coefficients can be estimated from the pseudo-BipoSH coefficients by 
inverting the above equation. 

Figure 3 shows BIPS for the foreground cleaned polarization maps recently 
released by WMAP filtered to retain power multipole between / ~ 5 to 15. The 
BiPS for CMB E-polarization maps at V and W from WMAP are compared 
to the mean and 1-a deviations of BiPS measurements on 100 simulated maps 
The appropriate galactic mask for polarization maps suggested by WMAP was 
employed for all maps. These preliminary results indicate strong violation of 
statistical isotropy in the polarization. However, the dependence on frequency 
suggests a non cosmological origin in foreground residuals in parts of the sky 
beyond galactic mask. Interestingly enough, c onsistent with conclusions of the 
WMAP polarization analysis, the BiPS analysis also indicates that the V band 
is 'cleaner' than the W band. 



4 Discussion 

The Bipolar representation of the statitical correlations underlying the CMB 
anisotropy and polarization allows a unprejudiced evaluation of statistical 



46c 47 ) whose elements 
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isotropy on largest observable cosmic scales. The null result of search for 
departure from statistical isotropy in the WMAP data provides strong evi- 
dence for cosmological principle and constrains non-trivial cosmic topology 
and ultra-large scale structure. It can potentially be used to constrain any 
other physics that violate global isotropy on cosmic scales. The null result has 
implications for the observation and data analysis techniques used to create 
the CMB anisotropy maps. Observational artifacts such as non-circular beam, 
inhomogeneous noise correlation, residual striping patterns, and residuals from 
foregrounds are potential sources of SI breakdown. 
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